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Towards the graviton from spinfoams: 
the 3d toy model 


Simone Speziale* 

Perimeter Institute, 31 Caroline St. N, Waterloo, N2L 2Y5, Ontario, Canada. 

Recently, a proposal has appeared for the extraction of the 2-point function of linearised quantum 
gravity, within the spinfoam formalism. This relies on the use of a boundary state, which introduces a 
semi“classical flat geometry on the boundary. In this paper, we investigate this proposal considering 
a toy model in the (Riemannian) 3d case, where the semi-classical limit is better understood. We 
show that in this limit the propagation kernel of the model is the one for the harmonic oscillator. This 
is at the origin of the expected l/£ behaviour of the 2-point function. Furthermore, we numerically 
study the short scales regime, where deviations from this behaviour occur. 


I. INTRODUCTION 

The spinfoam formalism Q is a candidate covariant approach to a non-perturbative quantisation of General Rela¬ 
tivity (GR). At present, it seems to provide a consistent background independent theory at the Planck scale, where 
it describes a (discrete) quantum geometry. However, the large scale behaviour is less understood. In particular, the 
formalism lacks a well-defined procedure to study the semi-classical limit, define particle scattering amplitudes and 
reproduce low-energy physics. Consider the pure gravity case: in the low-energy limit, we would expect to recover 
the conventional perturbative expansion described in terms of gravitons. The importance of obtaining this would be 
two-fold: on the one hand, it would provide fundamental ground for the emergence of General Relativity in the large 
scale limit of spinfoams, thus proving it to be a sensible theory of quantum gravity. On the other hand, it would 
open the way to the exploration of the corrections to the propagation of gravitons due to the microscopical quantum 
geometry, where new physics is potentially expected. 

There is a simple heuristic reason that shows the difficulty to recover conventional quantum field theory (QFT) 
from background independent approaches: the arguments of the n-point functions used in QFT are spacetime point 
coordinates, which are not defined in the absence of a background. Recently, a possible way out of this problem 
has been suggested in |2 • The key idea is to construct n-point functions by means of the propagation kernel, which 
provides an amplitude for the fields assigned on a boundary of spacetime. Since the value of the gravitational field 
provides the boundary with a metric, n-point functions can be defined with respect to this metric. In the deep 
quantum gravity regime, one is willing to consider arbitrary boundaries |^Q. However, to recover the behaviour 
of linearised quantum gravity, it seems as a hrst step sufficient to simply consider the usual setting of QFT, when 
the boundary is given by two spacelike hyperplanes, and the boundary values of the gravitational field are small 
perturbations over flat space, and which are awrmtotically vanishing. In this context, the propagation 
kernel W[h' ,h" ,T] for linearised GR was considered in @,0; explicitly derived in the temporal gauge in 0. 
The time T entering the kernel is the asymptotic proper time between the two hyperplanes. Indeed, it is a simple 
exercise to recover for instance the 2-point function Wabcd{x,y) by sandwiching W[h', h",T] with two one-particle 
states 'I'i[h']Qb(a;), ^i[h"]cd{y), and integrating over the boundary values and 

Wabcd{S,y;T) = 2 Jvh' Vh" W[h',h",T] 'i>o[h'] Ki,{0,x) 'I>o[h"] Kd{T,^, (1) 

where the normalisation constant N is given by the functional integral without the held insertions.^ Notice that both 
the kernel W and the vacuum state d'o (encoded in the kernel, see 0) are gauge invariant expressions; however, the 
held insertions and are not gauge invariant, thus the evaluation of the expression above requires an additional 
spatial gauge-hxing. Choosing for instance the Coulomb-like gauge dahab = 0, we have 

/ ^JTh _ 1 _ 2 

- 7 : - (OacObd + DadDbc - ^DabDcd\ ( 2 ) 
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^ Strictly speaking, this is only the advanced Green function A 4 -(ai,y). The Feynman propagator is obtained by taking 0(T)A+(a:,y) + 
e{-T)A-{x,y). 
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where n is the dimension of spacetime, and Dab = Sab — is the transverse projector. Different gauge choices result 
in different numerical factors in the tensorial structure. 

Following [2 , the 2-point function of linearised quantum gravity m can be translated into the spinfoam formalism 
using the expression 


Wf,„pa{x,y) = ^ W[s] «'o[s] {s\hp,a{x)\s) {s\hpa{y)\s). (3) 

S 

Here the sum is over spin networks s; the propagation kernel W[s] is provided by the explicit spinfoam model chosen, 
and it includes a sum over all spinfoams cr compatible with s. The boundary state 'I'o[s] is a weighted functional of 
spin networks peaked around those reproducing flat space in the semi-classical limit. 

This approach has been considered in Q for a proposal to extract the graviton from the (Riemannian) 4d Barrett- 
Crane model. The results there described are promising; indeed, considering a single component of the 2-point 
function, one does recover the correct dependence on the distance. However, there are many issues that still 
needs to be clarified. In particular, the role played by the Gaussian boundary state used to introduce the flat geometry 
on the boundary, and the way this geometry determines the distance dependence. 

Remarkably, these issues are present in a very similar way also in the 3d case. The 3d case has the advantage that 
the semi-classical limit is better understood, and there are no degenerate conhgurations plaguing the limit. In this 
paper, we propose to apply the strategy of [g in 3d, in the hope to clarify some of the issues. As it is well known, 
3d GR has no local degrees of freedom. Nonetheless, the 2-point function of the linearised quantum theory is a well 
defined quantity that can be evaluated once a gauge-fixing is chosen. For instance, it is given in the Goulomb gauge by 
m with n = 3; it shows the usual 1/p^ behaviour, which in 3d means a dependence 1/i over the spacetime distance. 
However, this quantity is a pure gauge, thus the quantum theory does not have a propagating graviton 

Notice also that in 3d the Newton constant G has inverse mass dimensions (in units c = 1); we define the 3d Planck 
length as ip = WirfiG. 


II. THE TOY MODEL 

The toy model we consider here is a tetrahedron with dynamics described by the Regge action, whose fundamental 
variables are the edge lengths ie- Since we are dealing with a single tetrahedron, all edges are boundary edges, and 
the action consists only of the boundary term, namely it coincides with the Hamilton function of the system: 

^R[4] = Y^^4 0e(4). (4) 

e 

Here the 9e are the dihedral angles of the tetrahedron, namely the angles between the outward normals to the triangles. 
They represent a discrete version of the extrinsic curvature M- Furthermore, they satisfy the non-trivial relation 

In this discrete setting, assigning the six edge lengths is equivalent to the assignment of the boundary gravitational 
field. In particular, let us consider the case when two opposite ed ges o f the tetrahedron have lengths a and 6, and the 
remaining four edges have all the same length, say c. Following |ll|. we know that the Hamilton function S[a,b,c] 
given by 0 defines a simple relativistic system, of which a, b and c are partial observables ^ 2 - That is, they include 
both the independent (time) variable, and the dependent (dynamical) variables, and the dynamics provides a relation 
amongst them. In the following, we choose c as the time variable. The classical dynamics of the system is described 
using the Hamilton function to read the evolution of b in a time c, given the initial value a. 

The quantum dynamics, on the other hand, is described by the Ponzano-Regge (PR) model ^3- the model, 
the lengths are promoted to operators whose spectrum is given by = ipC^{je), where the half-integer j 

labels SU(2) irreducible representations (irreps), and C'^(j) = j{j -|- 1) -I- | is the Casimir operator.^ In the model. 


^ Conventionally, the SU(2) Casimir is taken to be simply j(j + 1). However, it is defined up to an additive constant, and the value 1/4 
chosen here allows a matching with the original PR ansatz, ie = ^pOe + ^)- 
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each tetrahedron has an amplitude given by Wigner’s {6j} symbol for the recoupling theory of SU(2). Under the 
identifications a = ipC{ji), b = £pC{j 2 ), c = £pC{jo), the propagation kernel for our toy model is 

+ ( 6 ) 

e ^ 

It represents the amplitude for the (eigen)values of the boundary gravitational field. Notice that we can assume to 
have measured the quantity jo, and then interpret W[ji, j 2 , jo] as the transition amplitude between ji and j 2 in a 
time jo- 

The key property of the PR model which we will need in the following is that in the large j regime we have 

Emilia 

^Ihn IU[ji, J 2 , jo] ~ n + 1) + l) ’ 

where the argument of the cosine can be recognised to be {1/h times) the Regge action 0), with £e = £pC{je) (the 
constant factor j appearing in o clearly does not affect the classical dynamics). Usually, we would expect the 
semi-classical limit to be the exponential of the classical action. However, the presence of the two exponentials giving 
rise to the cosine can be understood as follows. The sign in front of the action is determined by the orientation 
of the manifold; it is reversed if we invert the orientation of M. However, the classical theory does not distinguish 
between the two orientation, thus the path integral sums over both. We will see below that only one contributes to 
the propagator. Therefore, the model admits a clear semi-classical limit, obtained by taking j i—> oo, £p i—> 0, such 
that the product £p j is finite. 


III. THE BOUNDARY STATE 

To make a link between this toy model and the 2-point function of GR, we consider two hyperplanes in 3d Euclidean 
spacetime, and we embed the tetrahedron as in FigUl Namely, we assume that we have measured the time interval 



FIG. 1: The dynamical tetrahedron as evolution between two hyperplanes. The labels give the physical lengths as a = £pC{ji), 
b = £pC{j2), c = £pC{jo), and T = £2 —t\ = cj-spi. 


T '■= £2 — ti, and we compute the correlator between perturbations of the length of the bottom edge, say £a, and 
perturbations of the top edge, £b. Because the two edges are embedded in fiat space, we can assume the coordinate 
vector £^ to be along the direction and £(( along the direction, as in FigUl With this setting, we are naturally 
led to consider the projection £'j{x)£^{x)£'^{y)£'l{y)Wij,^pa{x, y) of @, where x and y are the mid-points respectively 
of £a and £},'., therefore we restrict our attention to the following quantity. 


£l £l W^MT) ■■= £^Px)£:ix)£P{y)£Uy)W^,pPx, y) = £l £l fjp 


“p e 


(27r)2 2ujp 


( 8 ) 


where fx is a numerical factor depending on the choice of gauge (in fx = 2). 

In this “minimalistic” toy model, the boundary spin network s appearing in m is the boundary spin network of 
the tetrahedron. Since this is a tetrahedral graph with links in one-to-one correspondence with the edges of the 
tetrahedron, we can simply use jg for the spin network labels. However, we stress that we are considering a situation 
where we have measured the “time” variable T, thus the labels for the edges in the bulk are fixed to the value jof 


Consequently the conjugate variables, namely the bulk dihedral angles, are maximally spread. Notice that this is different from the 
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With a little trigonometry, we have T = c/-\/2- Therefore, the sum over s in 0 is truly simply a sum over ji and j 2 , 
and 0 is the propagation kernel to be used, 

ililWii 22 {x,y) = ^ W[jl,j 2 ,jo] 'i’Q[jl,j 2 ] {j 2 \e^ibhpa{y)\j 2 ). (9) 

ii J2 

The next ingredient entering 0 is the vacuum boundary state, j 2 ]- Its role is to peak the boundary spin 

networks around those reproducing flat spacetime. The natural choice is to take a Gaussian around a background 
value for the variables ji and j 2 - With respect to this, notice that any value of ji and j 2 is compatible with a flat 
tetrahedron. For simplicity, let us choose an equilateral tetrahedron: since the bulk edges are fixed to be jo by the 
time measurement, this means peaking the Gaussian around the value jo for both ji and j 2 - Furthermore, a key 
observation in Q is that the Gaussian should be peaked on the intrinsic as well as the extrinsic geometry of the 
boundary.'^ While the intrinsic classical geometry is represented by the edge lengths, the extrinsic geometry, on the 
other hand, must be fixed giving the boundary dihedral angles 0^- For the equilateral tetrahedron, the value of all 
dihedral angles is d = arccos(—|). Using jo and i?, we consider the tentative boundary state 

2 2 

J 2 ] = { “ f + ^)}- (10) 

^ 2—1 2—1 

Here a is the width of the Gaussian, which we leave as a free parameter for the moment. However, let us notice the 
following. In order for cni to represent a good vacuum boundary state, we require that in the semi-classical limit all 
relative uncertainties vanish, that is 


(Aj) 1 
{3) kVa 


m 

( 0 ) 



( 11 ) 


The semi-classical limit is obtained by taking jp 1 —> 00 . Since 1 ? is not affected by this limit, the conditions above 
imply a oc jV^ with 0 < r < 2. In spite of its simplicity, we show below that this Gaussian is enough to characterize 
the usual properties of the propagator over flat spacetime. 


IV. THE FIELD INSERTIONS 

The last ingredients entering © are the expectation values of the gravitational operator, such as (ji \^a^a^fj.v{x)\ji). 
A crucial point here concerns gauge invariance. Let us go back to O- As discussed in the introduction, the field 
insertions are not gauge invariant quantities, therefore an additional gauge-fixing is required, in order for o to be 
well defined. Precisely in the same way, a gauge-fixing has to be chosen to give a well defined prescription to identify 
(jil'^aCI^A»‘^(®)bi) ll^® P’11 model, and thus to evaluate 0. 

In particular, notice that we can choose in the continuum the gauge f-a^a^p.v{x) = 0; to be consistent with this 
choice, we have to set {ji\^a^a^p.v{x)\ji) = 0 in the PR model. Consequently, 0 is identically zero. The peculiarity 
of the 3d case is that this choice of gauge can be made simultaneously for all independent directions: the 3d graviton 
propagator is a pure gauge, and accordingly all the components can be set to zero. In the 4d case, on the other hand, 
the gauge choice = 0 cannot be made for all independent directions; in the end two components survive, 

the two degrees of freedom of a spin-2 massless particle. However, in the construction presented here we consider a 
single component of the graviton, thus we can not give a full proof of this peculiarity. To pursue this further, we 
need to extend the model to more tetrahedra, so that we can study the full tensorial structure. In fact, once we have 
the full tensorial structure, the counting of the number of physical degrees of freedom can be done using the Ward 
identities, which relate different diagrams through the symmetries of the theory. The relevant symmetry of 3d GR is 
the invariance under diffeomorphisms, which in the PR model is implemented by the so-called Pachner moves (see 
for instance ^3)’ these moves involve changing the number of tetrahedra in the triangulation. Therefore, in order to 


setting of |3, where all labels are allowed to fluctuate. However, a link between the two settings can be made by considering in |3 a 
Gaussian with different widths for the bulk and the top and bottom edges, and then studying the limit case when the former goes to 
zero. 

This is analogous to a Gaussian representing a coherent state for a free particle in QM, which is picked on both the position and 
momentum. 
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apply these moves and study the symmetries of the 2-point function, we have to extend the model to more tetrahedra, 
and we leave this issue open for further work. 

To study the large scale behaviour emerging from 0, we now proceed choosing a Coulomb-like gauge in which 
has components along the edge £q, so that its projections can be identified with non-trivial PR operators, such as 
the length operator introduced above. More precisely, since is a small perturbation around flat spacetime, 
we can write = (-a^aig^ivix) — 6fiu) ■= ipX^Xj^ — •^pC'^(jo), where ipC^{jo) is the background value of 

the length. Recall that X^X^ acts diagonally on an edge labeled by j, giving therefore, in this gauge, the 

semi-classical limit jo i—*■ oo gives 

(jilCC V(a^)lji) = ^l[C\h) - C^Oo)] ~ 2jo(ji - jo)^p. (12) 


V. GRAVITON PROPAGATOR 


We are now ready to perform the evaluation of 0. To shorten our notation, let us define Sji := ji — jo, i = 1,2. 
Using the explicit expressions 0, cni) and C2J, we have 


(I lUn22(r) = f E n + 1) m 


^T.i^Oi+i^T.iUi+i) r/y2 




(13) 


Jl ,32 e 


Before proceeding, let us stress the logic: we introduced a fiat boundary state, an equilateral tetrahedron with label 
jo; this enters into the choice of the Gaussian as well as into the definition of the field insertions. Then, to have a 
semi-classical behaviour of its geometry, we consider the large jo regime. In this regime, m is approximated by 

4jg Sji Sj2 

,/67rU(je) 


o*‘S'Rbe]+jf 


„-iSB.[3e]-ii 


„-§T,i^3i+i'^T,iUi + ^) 


(14) 


En(2je + 1) 

31,32 e 


This is the expression we now evaluate analytically to study the asymptotic behaviour of the 2-point function, which 
is to be compared with 0. As in Q , the Gaussian implies 5ji <C 1 in the sum, thus we can expand the Regge action 
around jo, 

-S'Rbe] = ~ E (jo + + E ^ 


Using 0, we have ^^\j^=jg = d. Introducing Gik = ^\j^=jg, we can write 


1 2 

S'Rb'e] ~ 4 d (jo -f -) -f E (ji 




(16) 


Gik can be calculated from elementary geometry; for the equilateral tetrahedron at hand, Gn = G 22 = 
Gi 2 = G 21 — If we insert this expansion into d, the two exponentials become 

^if+ii-&Uo + ^)+'2i'SJ2i (ti + 5 ) + 5 Ui.fc GikSjiSjk-i^Y,i ^3i _|_ g-*fUi.fc GikSjiSj^-^ Sj^ 


The first one is a rapidly oscillating term in ji, which vanishes when we perform the sum d, so that we can consider 
only the second term. Furthermore, notice that at this order in the expansion in Sji, we have U(je) ~ ^ijo), thus 
the volume does not enter the sum d; absorbing it in the normalisation, together with the phase e *1 4ii9(jo+2)^ 
we have 




1 

w 


4Jo4 E + 1) Sji6j2 

31,32 e 


(17) 


The sum can be easily performed by approximating it with a Gaussian integral. In fact, in the limit j 1 -^ 00 , ip 1 —> 0, 
we can write 


E + 1) ~ (2jo + 1)^ 

31,32 e 


(2ji + l)dji 


(2j2 + l)dj2 


(2jo)® 


djidj2, 
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and also the factor (2jo)® can be clearly absorbed in the normalisation. Changing variables from ji to Zi := 5ji, we 
have dji = dzi and 

el el vbii22(r) ~ 1 ajI 4 J dz^dz^ ziZ2 , (is) 

where Aij = aSij — iGij. At this point, we can fix the value of a to tune the correct vacuum boundary state. If we 
choose a = which is consistent with cu, we have 


Aij — ol 


i + ico%e 

- 51 ^ e + icote 

This can be easily recognised to be the kernel of an harmonic oscillator (see the appendix), with frequency 


(19) 


0^8 1 _2a 


( 20 ) 


Therefore, the vacuum boundary state together with the PR amplitude reproduce in the semi-classical limit the kernel 
of an harmonic oscillator; using the boundary geometry to evaluate 4^6 = (Jo^p)^, 


Wii22{T) ~ ^2 ^12^ 
Jo 


4 e*® 
io 2a 


8 

io^p 2a; 


( 21 ) 


The expression picks contributions from a single frequency. However, suppose now that our tetrahedron is part 
of a discretisation of spacetime. Then we expect to be able to sum over all frequencies admitted by the lattice, and 
becomes 


IPii22(T) 


8 

p 2a;; 


( 22 ) 


which coincides with (O, once we set fx = 8, and we make the usual lattice approximation f d'^p ~ (jolpy-^ Sp (where 
jo^p provides the characteristic length of the boundary geometry). 


VI. DEPENDENCE ON THE DISTANCE 


To study the dependence on the spacetime distance e of the 2-point function, it is convenient to consider its absolute 
value. From HU, we read that at large scales 


W{e) := |IFii22(^pjo)| = 


(23) 


thus reproducing the expected inverse power law of the 3d free graviton. 

However, we can now use the exact expression 0 to study the dependence on the distance at short scales. Notice 
that corrections to the inverse power law are expected to arise at short scales; in the context of conventional QFT, 
these are understood as graviton self-energies (see for instance ^3). In our approach, they are indeed present, and 
come from higher order terms in the expansion (trai as well as from departures from the asymptotic O- Here, we do 
not attempt an analytical evaluation of these corrections, but simply briefly report a numerical study: in Figl^we plot 
some numerical values of (US) when varying jg- We see that these are in very good agreement with the asymptotic 
behaviour m for jo > 50. For smaller values of jo, we begin to see deviations from the inverse power law. These 
deviations can be fitted by a series of inverse powers, which are the small scale quantum corrections. For instance, 
allowing powers down to —4, we obtain the following fitting function. 


w{e) 


+ ^2 0.74 


■4 


0.92 


(24) 


It should be added that these corrections are gauge-dependent quantities, just as the free propagator. In the 
conventional perturbative expansion of 3d GR, the self-energies vanish in the gauge in which the free propagator 
vanishes, so that the full propagator is a pure gauge. At the present stage of investigation, there is no reason to think 
that the spinfoam formalism should modify this picture, thus we expect these corrections to be again pure gauge 
quantities; namely, we expect them to vanish in the same gauge in which the free term vanishes. 












FIG. 2: The crosses are the numerical evaluations of cnj, plotted against the asymptotic behaviour ^ on a bi-logarithmic 
scale. The agreement is very good for jo > 50, while deviations appear for smaller values. 


VII. CONCLUSIONS 

We considered the recently appeared proposal for the extraction of the 2-point function of linearised quantum 
gravity from the spinfoam formalism. To clarify some of the geometrical issues, we focused on a 3d toy model where 
spacetime is discretised by a single tetrahedron. We showed that, upon properly fixing the width of the Gaussian used 
as the vacuum boundary state, the kernel for the space perturbations around flat space behaves as the kernel of an 
harmonic oscillator, with frequency given by where i = jo£p is the characteristic length of the boundary geometry. 
Therefore, the propagator can be reconstructed as a collection of harmonic oscillators, yielding the expected inverse 
power law in the large scale limit. Furthermore, we numerically studied the short scale behaviour, hinting that the 
1 j£ behaviour is corrected by a series of higher inverse powers. The proposal not only reproduces the expected large 
scale behaviour, but it is also useful to study the short scale quantum corrections. 

The main open issue is the extension of these results to many tetrahedra. On the one hand, this is an important 
check of the stability of the result; on the other hand, it would allow the study of the other components of the 
propagator, and thus recover its full tensorial structure. In particular, notice that the full tensorial structure is 
needed to prove in a more consistent way that the 2-point function computed here is just a pure gauge. 

We hope that the clarifications here presented will help the calculations in the 4d case. 
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APPENDIX A: 2-POINT FUNCTION FROM THE KERNEL FOR THE HARMONIC OSCILLATOR 

Consider an harmonic oscillator, with mass m = 1, in units fi = 1, and let us compute the 2-point function 
= (0|a:(ti)a;(t2)|0). This can be done in the canonical formalism, giving 


(Al) 



with t 2 — ti := T. However, it can also be computed starting from the propagation kernel VF[xi,X 2 ,T], as in |^. We 
have 


^ j dxidx2W[xi,X2,T] '^q[xi]xi'^o[x2]x2, 


(A2) 


where the normalisation is Af = f d:^ dx 2 a; 2 , T] 'I'o[®i] ^o[a; 2 ], and 'I'oN is the vacuum state. Using the well 


rt 

known expressions (see for instance [ij) 


^oN = W[xuX2,T] = 

V TT 


UJ {xf+x^)cosujT-2x-^X2 

-g * 2 sinwT 


27rz sina;T 


(A3) 


we obtain 


G{ti,t2) = ^ J dxidx2XiX2e ^ 


(A4) 


with 


1 + i cot uT — 


^ij — ^ 


sin uiT 


sin ljT 

1+1 cot ioT 


(AS) 


The expression (IA4ll is a Gaussian integral that can be readily evaluated, giving 


G{ti,t2) = Ai2 = 


iujT 


(A6) 


Up to the vacuum energy contribution the result coincides with the canonical evaluation m- 
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